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We present measurements of second- and higher-order intensity correlation func- 



tions (so-called Hanbury Brown and Twiss experiment) performed at the free- 
electron laser (FEL) FLASH in the non-linear regime of its operation. We demon- 
strate the high transverse coherence properties of the FEL beam with a degree of 
transverse coherence of about 80% and degeneracy parameter of the order 10 9 that 
makes it similar to laser sources. Intensity correlation measurements in spatial and 
frequency domain gave an estimate of the FEL average pulse duration of 50 fs. Our 

> 

measurements of the higher-order correlation functions indicate that FEL radiation 

in ; 

obeys Gaussian statistics, which is characteristic to chaotic sources. 
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Hanbury Brown and Twiss in their pioneering experiments demonstrated that one 
can get fundamental information on the statistics of light sources by measuring intensity 
correlations at two separated spatial positions. Originally designed as a robust method to 
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determine the size of stars, these experiments initiated developments in the field of quan 



turn optics [3(. Statistical properties of thermal sources 



microcavities [5], and, recently, Bose-Einstein condensates 



lasers 



4J, semiconductor 



3] have 



i been studied using 



this technique. The recent advent of x-ray free-electron lasers (FELs) 



11] with their un- 



precedented peak brilliance and ultrashort pulse duration has opened the route to a number 



of spectacular ground-breaking experiments including femtosecond nanocrystallography 12 1 
and single particle coherent imaging [13I . 14j . Many of these experiments exploit the high 
degree of coherence of the FELs. 



Coherence is the defining 



within statistical optics 15 



eature of a laser source and is described by correlation functions 



16]. The first-order correlation properties of FEL sources have 



been extensively investigated recently 17H22]. It was experimentally demonstrated that 
FELs based on the self-amplified spontaneous emission (SASE) process have a high degree 
of transverse coherence but poor temporal coherence. To get a more detailed picture of the 
statistical properties of these sources, higher-order field correlations must be studied. These 
can be explored, for example, by utilizing intensity correlation measurements in a Hanbury 
Brown and Twiss (HBT) experiment. In this letter we present measurements of second- and 
higher-order intensity correlation functions at the free-electron laser FLASH 8]. 

The core idea of the HBT experiment [l, 2] is to determine the normalized second-order 
intensity correlation function 

(/(ri) ■ I(r 2 )> 



<? (2) (ri,r 2 ) 



(1) 



<J(ri))</(r 2 )>' 

by measuring the coincident response of two detectors at separated positions and r 2 
(see for review [23]). In Eq. (CQ), I(ti), I(r 2 ) are the intensities of the wavefield, and the 
averaging is done over a large ensemble of different realizations of the wavefield. It is well 
established that chaotic light can be described in the frame of Gaussian statistics 15] and 
is completely determined by the first-order correlation function in spatial domain known as 
the normalized spectral degree of coherence (SDC) r 2 ). It is defined as [15] /i(ri, r 2 ) = 
W{ri, r 2 )/ ^/S(r 1 )S(r 2 ), where H / (r 1 ,r 2 ) is the cross spectral density function and S(r) is 
the spectral density. 

The intensity correlation function then reduces to (see the Appendix A for details) 



g^(r 1 ,r 2 ) = l + ( 2 (D UJ )\ f i(r 1 ,r 2 ) 



(2) 



3 



where Cz(Du) is the contrast function that strongly depends on the bandwidth D u of the 
radiation. It was earlier demonstrated [24J that for stationary chaotic sources the contrast 
Ci{Dj) is determined by the ratio t c /T, where r c = 2ir /D u is the coherence time of the 
wavefield and T is the time resolution of the detectors. Equivalently, the contrast ^iD^) 
determines the number of longitudinal modes Mt = l/^iDj) for a chaotic source. By 
definition the SDC |/x(r 1? r 2 ) | < 1 and intensity correlation function g^(ri,r 2 ) < 2. 

For the full statistical description of the wavefield the n-th order correlation functions 
can be introduced 



15] 



9 (ri, — r " ) - re, </(*»• (3) 

For chaotic light q^fri , . . . , iv) is completely described by the SDC /x(ri,r 2 ) due to the 

15] . A comparison between correlation functions of different 



Gaussian moment theorem 



orders determines whether the field obeys Gaussian statistics or not. In particular, for 
chaotic light the n-th order correlation function g^(r,r, . . . , r) is equal to n! for a single 



longitudinal mode 



15|- 



For pulsed sources the intensity correlation measurements are naturally gated by the 



pulse duration T 



241 ] . FEL sources, with pulses of few tens of femtoseconds, are ideally 



suited for intensity correlation measurements. According to FEL theory 25| these sources 
should obey Gaussian statistics in the linear and deep non-linear regime of operation. In 
these conditions the correlation function g^(ri,v 2 ) has the form of Eq. (T5]) and provides 
access to the transverse coherence properties of an FEL as well as to its pulse duration 24]. 

The experiment was carried out at FLASH that was operated with six undulator modules 
and a total undulator length of 30 m. The electron bunch charge was 600 pC, and the electron 
energy 1.08 GeV resulting in a photon wavelength of A=5.5 nm. The average photon pulse 
energy was about 110 /xJ, which corresponds to about 3-10 12 photons per pulse at this photon 



energy. The measurements were performed at the PG2 beam line [26j, |27j (see Fig. [I]). The 
optical system focuses the beam at a distance of 71.5 m downstream from the undulator exit. 
A monochromator comprised of a plane grating (PG), collimating (Ml) and focusing (M2) 
mirrors, and an exit slit with variable slit width was utilized to modify the bandwidth. The 
plane grating has a line density of 200 lines/mm and was tuned to its third order resulting 
in a dispersion in the exit slit plane of 0.64 eV/mm. The resolution of the monocromator 
was E=7 meV, which corresponds to an energy bandwidth of AE/E = 3.1 ■ 10~ 5 . Exit slits 
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FIG. 1: (Color online) (a) Scheme of the experiment. FEL radiation is generated in the undulator 
and is transmitted through the beamline including three mirrors (Ml, M2, M3), the plane grating 
(PG), and exit slit. Intensity profiles of individual femtosecond pulses are measured at the detector. 
(b,c) Typical single pulse intensity profiles and (d) an average over 2 • 10 4 pulses for a bandwidth 
of AE/E = 0.8 • 10 -4 . (e-f) Projections of the pulse intensities along the vertical direction. The 
scale bar is 0.5 mm long. 

of 30 /im, 60 fim, 150 /im, 300 /im and 500 /im in size were used to select the spectral width. 
These values correspond to energy bandwidths AE/E of 0.8 • 10~ 4 , 1.7 • 10 -4 , 4.0 • 10~ 4 , 
0.8 ■ 10~ 3 , and 1.4 ■ 10~ 3 . 

In contrast to the original HBT experiment [l, 2| with two separated detectors, we utilize a 
pixel detector for correlation measurements. It allows to determine second- and higher-order 
correlation functions at all separation points simultaneously. An in-vacuum CCD (Andor 
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FIG. 2: (Color online) Intensity correlation analysis. (a,b) Intensity correlation function g( 2 \xi, x%) 
for a bandwidth of 0.8 • 10 -4 (a) and 1.4 • 10~ 3 (b). (c,d) Intensity correlation function g^(Ax) 
taken along the white line in (a,b) around the center of the beam. Insets in (c,d) show intensity 
fluctuations g( 2 \x,x) = (l 2 (x)} / (I(x)} 2 taken along the blue line in (a,b). The error bars in (c,d) 
are obtained by statistical analysis of 20 individual sub-ensembles (10 3 shots each) from the whole 
set of 2 • 10 4 shots. The number of points displayed in this figure with error bars is reduced for 
better visibility. 

Ikon, 2048 x 2048 pixels, each 13.5 /im x 13.5 /im in size) was positioned at a distance 
of 3.3 m behind the focus of the beamline. The detector was operated at a repetition rate 
of 10 Hz with a region of interest 2047(H) x 460 (V) pixels and binning by five pixels in 
the vertical direction. A silicon nitride film 14 mm x 14 mm in size and 1 /im thick was 
positioned approximately 30 cm upstream from the camera to attenuate the beam. The 
transmission of the film at this energy was 4 • 10~ 6 . About 2 • 10 4 intensity profiles were 
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FIG. 3: (Color online) Statistical properties of FEL light as a function of the bandwidth, (a) 
Contrast (points) C2(A>) = fi f(2) (°) ~ 1 

as a function the coherence time r c . The error bars are 
same as in Fig. 2 (c,d). Theoretical fit (solid line) correspond to a pulse duration of T = 47±21 fs. 
(b) Measured single pulse (blue line) and average (black line) spectra. Gaussian fit (red line) gives 
an FEL bandwidth of AE/E = 6.7 • 1CT 3 FWHM. The width of the largest exit slit used in the 
experiment and corresponding to the bandwidth AE/E = 1.4- 10~ 3 is shown by the gray region, (c) 
Intensity correlation function </ 2 )(wi, 0J2) in the spectral domain, (d) Intensity correlation function 
g( 2 )(Au;) (red circles) taken along the white line in (c). Gaussian fit (blue line) gives an average 
pulse duration of 27 fs. 

recorded for each monochromator setting. Dark images were subtracted from the measured 
intensity distribution, and occasionally occurring negative values were set to zero. Typical 
single pulse and averaged intensities are shown in Fig. [T] (b-d). Shot to shot fluctuations 
in these intensity profiles, as a consequence of the SASE process are clearly visible. The 
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FIG. 4: (Color online) Higher order correlation functions for a bandwidth of 0.8 • 10~ 4 (a, b) 
and 1.4 • 10 -3 (c, d). (a,c) Third-order correlation function g( 3 '(x\, X2) (shown in the insets) and 
g( s )(Ax) (black solid line) taken along the white line in the insets. (b,d) Fourth-order correlation 
function g^ A \x\, X2) (shown in the insets) and g^(Ax) (black solid line) taken along the white line 
in the insets. The same functions obtained under assumption of the Gaussian statistics are shown 
as blue dashed lines. The number of points displayed in this figure with error bars is reduced for 
better visibility. Error bars are evaluated in the same way as in Fig. [2j 

intensity profiles were averaged along the vertical direction, which is the dispersive direction 
of the monochromator (see Fig. [1] (e-g)). The intensity correlation analysis was performed 
in the horizontal direction. 

The normalized second-order correlation function g( 2 \xi, X2) for a narrow bandwidth of 
0.8 ■ 10~ 4 is shown in Fig. [2] (a,c). Remarkably, it reaches the maximum value of two at small 
separations (Fig. [2] (c)), which indicates that the contrast ^{D^) at that monochromator 
setting is close to one. This contrast is significantly higher than at synchrotron sources 28- 
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30J, where it did not exceed 0.3. A Gaussian fit, exp ((— Ax 2 )/(2Z 2 )), to the second-order 
correlation function g( 2 \Ax) as a function of the separation Ax = x 2 — X\ around the center 
of the beam (see Fig. [2] (c)) provided a transverse coherence length of l c = 0.93 ± 0.04 mm. 
This value is substantially larger than the measured beam size in horizontal direction (0.45 
mm (FWHM)), indicating a high coherence of the beam. We quantified it by evaluating the 



degree of spatial coherence (s defined as [3l|, |32| (s = f \W(xi, X2)\ 2 dxidx 2 /{J S(x)dx) 2 and 
obtained — 0.78 ± 0.01, which is in a good agreement with the Young's double pinhole 
measurements at FLASH 



211 ] . Importantly, we can estimate the degeneracy parameter 



151, 



which is the number of photons in a single mode. Our estimates, obtained by integrating 
the total flux on the detector at the narrow bandwidth, yield a value of 10 9 comparable to 
Ref. j^lj], which is significantly higher than at any synchrotron sources. 

The second-order correlation function g^(xi, x 2 ) for a larger bandwidth of 1.4 • 10~ 3 
is shown in Fig. |2]^b,d). This larger bandwidth is equivalent to a shorter coherence time, 
and the contrast of g^(xi,x 2 ) is reduced as expected from Eq. (|2]). Unexpectedly, at 
these conditions we observed an oscillatory behavior of the correlation function g( 2 \Ax) 
(see Fig. EJ^d)). This may originate from the contribution of two independent sources in 



the lasing conditions of FLASH as discussed in Ref. [33] . We want to remark that such 
fine features of FEL lasing would be difficult to observe by other means. For example, our 
averaged spectral measurements (see Fig. [3]^b)) do not show any notable features. 

Results of our experiments also indicate that the statistical properties of the FEL beam 
are not spatially uniform. This is well seen by the inspection of the beam fluctuations 
g( 2 \x, x) = (I 2 (x)) I (I{x)) 2 along the beam profile (see insets in Fig. |2]^c,d)). They vary 
significantly both for small and large bandwidths. 

We also analyzed the contrast in the center of the beam ^(D^) = g^ 2 '(0) — 1 as a function 
of the coherence time r c (see Fig. [3]^a)). It has a linear dependence as a function of r c at a 
large bandwidth and reaches saturation at a small bandwidth, as described by Eq. (|2]). A 
theoretical fit to the contrast values yields an average pulse duration of 47 ± 21 fs (FWHM) 
(see the Appendix B for details) 

To get an independent estimate of the pulse duration, we measured the second-order 
correlation function g^(oJi, co 2 ) in the spectral domain as a function of two frequencies (see 



Fig. [3]^c)) 34j, [35(]. About 1.5 ■ 10 4 single shot spectra were recorded after the intensity cor- 



relation measurements with a detector at the position of the exit slit of the monochromator. 
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This detector is comprised of a scintillating screen (YAG:Ce 0.2%) and an intensified CCD 
(Andor iStar, DH740), equipped with a lens. The effective pixel size of the detector in the 
exit slit plane was 19.4 /im with the point spread function estimated to be about two pixels 
(FWHM). The detector was operated at a repetition rate of 10 Hz. Line profiles, obtained 
from the central part of the beam in the horizontal direction corresponding to about 10 
% of the beam FWHM were analyzed. These line profiles were background corrected by 
subtracting a constant offset, which was obtained from the averaged spectrum. The oc- 
casionally occurring negative values in these line profiles were set to zero. A Gaussian fit 
exp (—Au> 2 ■ T 2 /2) was used to determine the pulse duration T from the measurements of 
the second-order correlation function g( 2 \A<jj) (see Fig. EJ^d)). According to this analysis we 
obtained an average pulse duration of 27 fs (FWHM), which lies in the uncertainty range of 
the previous measurements. In general, our estimates of the pulse duration can be shorter 
than the intrinsic pulse duration in time domain due to a possible frequency chirp of FLASH 



pulses [35 



To get an insight into the photon statistics of FEL pulses we studied intensity correlation 
functions of order higher than two. Higher-order correlation functions g^(xi,X2,Xs) (with 
x% = 0) and g^(xi,x 2 , x 3 , x 4 ) (with x 3 = xi/2 and x 4 = x 2 /2) are presented in the insets of 
Fig. H]for two different bandwidths as functions of X\ and x%. In Fig. H]the same quantities 
are presented as a function of Ax around the center of the beam with x\ = —Ax/2, x 2 = 
Ax/2. The higher-order correlation functions were also calculated under the assumption 
of Gaussian statistics (see the Appendix A for details) from the measured second-order 
correlation function g <y2 \x\,X2) (see Fig. H|). The excellent agreement between these two 
curves confirms the chaotic nature of the FLASH source at these operation conditions. In 
particular at a small bandwidth we observed a factorial behavior of higher order correlation 
functions at zero point separation g^ n '(r, r, . . . , r) ~ n\ (see Fig. H^a,b)), which is typical 
for Gaussian statistics. At a wider bandwidth the maximum value is substantially reduced 
compared to a narrow bandwidth case (see Fig. H|c,d)). 

In summary, we have presented intensity correlation measurements at FLASH which 
provide a simple, robust and versatile tool for monitoring basic beam properties of FELs 
including the degree of spatial coherence, average pulse duration, and details of the photon 
statistics. We obtained an averaged pulse duration of 50 fs, degree of transverse coherence 
of about 80%, and degeneracy parameter of the order of 10 9 . Such values are similar to laser 



10 



sources and were never observed at these wavelengths at conventional synchrotron radiation 
sources. However, our measurements of the higher-order correlation functions indicate that 
present FEL sources based on the SASE principle are essentially chaotic sources obeying 
Gaussian statistics. 

An interesting further application of the methods developed in this work would be the 



study of seeded FEL sources [37|, [38(. An intriguing question is whether seeded FELs are 



fully coherent sources in all orders according to Glauber [3] and in this way are equivalent 
to conventional single mode lasers, or if they obey Gaussian statistics like SASE FELs. 
We could also foresee that intensity correlation analysis might be applied to study the 
dynamics of ultrafast processes at FELs. For example, the intensity correlation analysis of 
the Coulomb explosion in single molecule imaging experiments [39| could provide detailed 
information about the disintegration of these molecules on a femtosecond time scale. 
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Appendix A: Higher order intensity correlation functions 



The normalized intensity correlation function is defined as 

(2) _ , (-f(r0/(r 2 )) 



^•*> = #TOy- (A1) 
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where r 1; r 2 are positions in space. In our implementation of the Hanbury Brown and Twiss 
experiment the intensity is given by 

oo 

J(r)= J \T(r,co)\ 2 \E(r,u)\ 2 du, (A2) 

— oo 

where E(r,u) is the radiation field in the space-frequency domain and T(r,u) is the trans- 
mission function of the monochromator. The average (• • • ) in Eq. (lAlj) is performed over 
an ensemble of different FEL pulses. In our experiment the intensity correlation function 
<?( 2 )(ri,r2) represents a coincidence measurement between different pixels of the detector. 

Below we consider only the horizontal direction x and a monochromator with a trans- 
mission function T(u) that does not depend on the position. Substituting Eq. (IA2j) into 
Eq. (IAip and interchanging the order of average and integration we find 



/ / \T(u 1 )\ 2 \T(u 2 )\ 2 (E^x 1 ,u 1 )E%x 2 ,co 2 )E(x 2 ,co 2 )E(x 1 ,co 1 ))du 1 du 2 



g {2) (x 1 ,x 2 ) 



f f ir^OHTMPdE^^oi^-d^^^^i^d^d^ 

— oo — oo 

(A3) 

For fields obeying Gaussian statistics, higher-order correlations of the field can be expressed 



through the first-order correlation functions using the Gaussian moment theorem 15] 

(E*(xi,ui) ■ ■ ■ E*(x n ,u n )E(x n ,u n ) ■ ■■E(x 1 ,ui)) 

= y^ j (E*(xi,uj 1 )E*(x^ 1 ),u n ( 1 ))) ■ ■ ■ (E*(x n ,u) n )E*(x n ( n ),u} n ( n ))). 

Here 7r denotes a summation over all n\ possible permutations of the subscripts. Applying 
this theorem to Eq. (IA3I) with n = 2 yields 



/ / iT^mu^lW^ixux^u^u^d^dco, 

g^(x 1 ,x 2 ) =1 + , (A5) 

/ / \T(u 1 )\ 2 \T(u 2 )\ 2 S(x 1 ,u 1 )S(x 2 ,uj 2 )duj 1 du 2 

— oo — oo 

where 

W^{xt, x 2 , u u tu 2 ) = (E*{x 1 ,co 1 )E{x 2 , co 2 )) (A6) 

is the first-order correlation function in space-frequency domain, or cross spectral density 
Q and 

S(x,u) = W w (x,x j uj,uj) (A7) 
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is the spectral density. 

With the assumption that the cross spectral density and spectral density can be separated 
into its spatial and spectral part 

W (1) (X 1 ,X 2 ,UJ 1 ,UJ2) = W(X 1 ,X 2 )W(U 1 ,UJ2), 

(A8) 

S(x, u) = S(x)S(oj) 
the intensity correlation function reduces to the form 

g {2 \xi, x 2 ) =1 + ( 2 {D U ) ■ x 2 ) | 2 , (A9) 

where 

W(x 1 ,x 2 ) . 

Ufa, x 2 ) = (A10) 

y/S(Xi)S{X2) 

is the normalized spectral degree of coherence. In Eq. (1A9h the contrast £ 2 (D U ) of the 
second-order correlation function is defined as 

oo oo 

C2p«) = — ~ 2 (AH) 

(/ \T(u)\ 2 S{u)duj) 

— oo 

and D w is the spectral width of the monochromator. The functions fi{uJi,uj 2 ), and S(oj) are 
defined accordingly from W{u)\,lo 2 ). 

The third-order correlation function is given by 

Substituting Eq. (1A2|) into Eq. (IA12I) . applying Gaussian moment theorem (1A4[) . and using 
approximation flA80 we find 

g i3 \xi,x 2 ,x 3 ) =1 + ( 2 {D U ) ■ (\fj,(xi,x 2 )\ 2 + \n(x 2 ,x 3 )\ 2 + \fi(x 3 , x^f^j 

+2( 3 (D U ) ■ Re(jj l (x 1 ,x 2 )fj l (x 2 ,x 3 )fj l (x 3 ,x 1 ) 

where ( 2 (Du) has been defined in (1A11I) and the contrast ( 3 (Du) for the three-point corre- 
lation function is given by 

oo oo oo 

/ / / |T( Wl )| 2 |T( W2 )| 2 |T( W3 )|W( Wl , W2 )W( W2 , W 3)^( W3 ,Wi)da; 1 da; 2 du;3 



(A13) 



-oo — oo — oo 



/ \T{u)\ 2 S{u)du ' 



(A14) 
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In a similar way the fourth-order correlation function 

(4)/ \ (I{x 1 )I(x 2 )I{x 3 )I{x i )) 

g { ) {x 1 ,x 2 ,x 3 ,x 4 ) = j- A15 
{I{x 1 ))(I{x 2 ))(I{x 3 ))(I{x^) 

can be expressed through the first-order correlation functions in the frame of Gaussian 
statistics 

g (4) (x 1 ,x 2 ,x 3 ,x 4 ) = 1 + (2(D U ) ■ (\^{xi,x 2 )\ 2 + |/i(xi,x 3 )| 2 + |/i(xi,x 4 )| 2 

+ Ci( D u) ■ (\^2,X 3 )\ 2 + | / u(x 2 ,X 4 )| 2 + |/i(x 3 ,X 4 )| 2 j 

+ Ci( D u) 2 • ^|//(xi,x 2 )| 2 |/i(x 3 ,a;4)| 2 + |/i(xi,x 3 )| 2 |/i(x2,x 4 )| 2 + \n(x 1} x 4 )| 2 |/i(x 2 , x 3 )| 2 

+ 2(3(1^) • Re (/i(x 2 , x 3 )/j,(x 3 , x 4 )fi(x 4 , x 2 ) + //(xi, x 3 )fi(x 3 , x 4 )/i(x 4 , x x ) 

+ 2C 3 (£> W ) • Re^/i(xi,x 2 )/i(x 2 ,x 4 )/i(x 4 ,xi) + fi(x 1 ,x 2 )fi(x 2 ,x 3 )fi(x 3 ,x 1 ) 

+ 2C 4 (L> W ) • Re^(xi,x 2 ) / u(x 2 ,x 3 ) / u(x 3 ,x 4 ) / u(a; 4 ,a:i 

+ 2C 4 (£) W ) • Re^/i(xi,x 3 )/i(x 3 ,x 2 )/i(x 2 ,x 4 )/i(x 4 ,xi 

+ 2C 4 (£> W ) • Re^(xi,x 3 ) / u(x 3 ,x 4 ) / u(x 4 ,x 2 ) / u(x 2 ,xi 

(A16) 

where ^(D^) and ( 3 (D U ) have the same meaning as in Eqs. (1A11I) and ()A14j) . The contrast 
C4(A*>) for the four-point correlation function is defined as 

f 1 n \ _ / f( u i> w 2)/(w2, w 3 )/(w 3 , W4)/(w 4 , u 1 )duj 1 du 2 duj 3 du 4 

(/|T( W )|^( W )(L;) 

where = T*(wj)T(wj)W / (co'i, Wj). The quantities (2^), CsC^^) , and C4.(Aj) are rea l 

numbers. The phase of the cyclic product of SDC in Eqs. (1A13I) and ( 1A16I) vanishes in most 
practical cases, and was thus neglected here. This includes fully coherent and incoherent 
beams, and Gaussian Schell- model beams [l5 |. 

To investigate whether FEL radiation obeys Gaussian statistics, we determined the mod- 
ulus of the spectral degree of coherence \fi(xi,x 2 )\ and contrast C2(Aj) from the second- 
order correlation function measurements according to Eq. (1A9I) . We then compared the 
third-order correlation function determined from the experiment (see Eq. ( 1A12I) ) with the 
expression ( 1A13I) . The same was done for the fourth-order correlation function determined 
from the experiment using Eq. (1A15j) and compared with Eq. (1A16j) . The parameters ^(D^) 



and Ca{D u ) in Eqs. (IA13I) and (|A16[) were calculated as a function of ^(Dj) assuming a 
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Gaussian form of the monochromator transmission function T(u) = exp (—u 2 /[A ■ D%\) and 



Gaussian Schell-model pulses 

Cs(A0 



40j |. In the frame of this model we obtained 



3 + 



and C^A. 



(A18) 



Appendix B: Determination of the pulse duration from the intensity correlation 

measurements 

To determine the pulse duration we analyzed the contrast (,2{DJ) obtained from the 
intensity correlation measurements performed with different monochromator settings. The 
transmission function of the monochromator was considered in the following form 



T(u) 



(Bl) 



1 for \u\ < 
for |w| > DJ2 

The FEL pulses before the monochromator were taken in the form of Gaussian Schell-model 



pulses [40j. In this model the spectral part of the cross spectral density V4 7 (wi,w 2 ) (see 
Eq. fTASjl ) is given by [40] 

(ui - UJ ) 2 + (u) 2 - UJ ) 2 (u 2 -uji) 21 



W(ui, 0J2) = exp 



AO? 



2tt 2 



(B2) 



where Q is the spectral width of the pulse before the monochromator, Q c is the correlation 
width of the spectrum, and uq is the mean angular frequency. These parameters can be 
expressed through the r.m.s. values of the pulse duration T and the coherence time T c of 



the pulse before the monochromator 



40J 



ft 2 = ^7 + 



rp2 



and Q c = -^Q. 



(B3) 



In the conditions of our experiment the coherence time T c before the monochromator is 
much smaller than the pulse duration T. Then the correlation width, Q c , in the spectral 
domain is well approximated by 



(B4) 



Under the assumption that the transmitted bandwidth is much narrower than the bandwidth 
of the incoming FEL radiation Q, (see Fig. 3(b) of the main text) the contrast C2(Aj) 
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as a function of the transmitted bandwidth can be calculated by integrating Eq. (lAllj) 
analytically 

Ca(A,) = ^erf {D W T) + ^-L^ {exp [-(A.T) 2 ] - 1} . (B5) 

According to Eq. (IB5I) the contrast C,2{D U ) approaches the value of one for D^T <C 1 and it 
has a linear behavior as a function of 1/D U when condition D^T ^> 1 is satisfied. According 
to the definition of the coherence time r c = we obtain in this limit Cs(Aj) ~ r c/T. 

The pulse duration was determined by fitting Eq. (1B5j) to the contrast ^(D^) for different 
measurements, corresponding to different bandwidths D u of the monochromator settings (see 
Fig. 3(a) of the main text). 



Appendix C: Results of the measurements from all monochromator settings 

The results from all measured monochromator settings are summarized and presented in 
Table 1. 

Table 1 



AE/E 


Beam FWHM 


Coherence Length 


De^ 


;ree of coherence 


Contrast 


Number of modes 


0.8 • 10~ 4 


0.45 mm 


0.93 ± 0.04 mm 




0.78 ±0.01 


1.03 ±0.05 


1.5 ±0.1 


1.7- 10~ 4 


0.45 mm 


0.93 ± 0.03 mm 




0.77 ±0.01 


0.97 ±0.06 


1.7 ±0.1 


4.0 • 10~ 4 


0.46 mm 


1.13 ±0.06 mm 




0.89 ±0.01 


0.76 ±0.04 


1.4 ±0.1 


0.8 • 10~ 3 


0.46 mm 


0.86 ± 0.03 mm 




0.82 ±0.01 


0.31 ±0.01 


3.5 ±0.2 


1.4- 10~ 3 


0.46 mm 


0.53 ± 0.04 mm 




0.72 ±0.01 


0.14 ±0.01 


7.8 ±0.4 



We also present the second-order correlation function for all measured monochromator 
settings in Fig. [5j The histogram distributions of the total intensity per pulse are also 
shown in the same Figure. Experimental results for the third and fourth-order intensity 
correlation functions as well as comparison with the expressions ( 1A13j) and (1A16j) for different 
monochromator settings are presented in Fig. [6l 
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FIG. 5: Left and middle column: Same as in Fig.[2]of the main text for all measured monochromator 
settings. The monochromator bandwidth AE/E was (a) 0.8 • 10 -4 , (b) 1.7 • 10~ 4 , (c) 4 • 10 -4 , (d) 
8 • 10~ 3 , and (e) 1.4 • 10~ 3 . Right column: Histograms of the total intensity per pulse (red bars). 



Theoretical fits with the gamma function 
(spatial and temporal) are also shown. 



311 ] (black lines) and the number of total modes M 
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FIG. 6: Same as in Fig. 4 of the main text for all measured monochromator settings. The monochro- 
mator bandwidth AE/E was (a) 0.8- 1(T 4 , (b) 1.7- 1(T 4 , (c) 4-10" 4 , (d) 8-l(T 3 , and (e) 1.4- 10~ 3 . 
Third-order correlation functions (left column) fourth-order correlation functions (right column). 
In the insets two-dimensional distributions of the corresponding correlation functions are presented. 
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